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Abstract 

The (l+l)-dimensional bosonization relations for fermionic mass terms are derived 
by choosing a specific gauge in an enlarged gauge- invariant theory containing both 
fermionic and bosonic fields. The fermionic part of the generating functional subject 
to the gauge constraint can be cast into the form of a strongly coupled Schwinger 
model, which can be solved exactly. The resulting bosonic theory coupled to the scalar 
sources then exhibits directly the bosonic counterparts of the fermionic densities ^/S^ 
and V'- 
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In a previous paper [Q], we have shown how the Abehan bosonization relations for the 
vector and axial vector currents in two dimensions @, ^ Q can be derived from a novel 
perspective, based on a new local "bosonization gauge symmetry" 0. The idea is, briefly 
stated, that bosonic and fermionic formulations of the same theory should be understood 
as two different gauge fixings of a larger gauge-invariant action containing both bosonic and 
fermionic fields. Equivalence between certain purely bosonic and purely fermionic formula- 
tions should then amount to the usual gauge-fixing independence of S-matrix elements in 
gauge theories. Although this idea is appealing, it is a non-trivial matter to show in detail 
that it can be carried through. One needs two ingredients: (a) The local bosonization gauge 
symmetry and the "larger" gauge-invariant action, and, (b) a smooth interpolating gauge- 
fixing function which can bring one continuously from a "boson gauge" to a "fermion gauge" 
by tuning the gauge-fixing parameter. 

In the usual path- integral approach to bosonization P, |^, |^ , it is clear that chiral fermion 
determinants play a fundamental role. The idea of [|T| was therefore to start with a fermionic 
theory, and then promote chiral rotations to a local chiral gauge symmetry. Using the general 
scheme of f^, this can be done without changing the physical content of the theory through 
the introduction of a collective field 9{x). This collective field - in essence the chiral phase 
of the fermions - turns out to be the bosonized field in a suitable gauge. In we showed 
how to introduce a particular gauge- fixing function depending on a gauge parameter 
A in such a way that one interpolates smoothly between purely fermionic formulations 
(A = 0) and purely bosonic formulations (A = 1). By comparing the couplings to external 
sources we recovered the usual bosonization relations for the currents. One of the most 
interesting aspects of this gauge-invariant approach to (l + l)-dimensional bosonization is 
that it demonstrates that these known equivalences are but two extreme cases of a continuum 
of equivalent theories that contain, in general, both fermions and bosons interactively. For 
this reason we dubbed our scheme smooth bosonization. 

Although the current bosonization relations ip'^^ip ~ and ip'jfj.75'4' ~ '^~^^'^d^9 

can be established in this manner by a fairly direct route, the derivation of the much more 
subtle mass-term bosonization relations 

i/jt/j ^ M cos{2^/7^9) and ip^y^t/j M sm{2^/^T9) (1) 

was only outlined in The purpose of this letter is to fill in this gap. 

Since we have already shown in [|I|] how to treat external vector and axial vector current 
sources, we shall only consider the generating functional 

= ^p{x){i0 + M+{x)P+ + M^{x)P^)ip{x), (2) 

which for constant sources M± just corresponds to having scalar and pseudoscalar mass 
terms. Here P± = (1 ± 75) are the usual chiral projectors. 

Although global chiral symmetry is broken explicitly by the sources, we still have what 
at first sight looks like a trivial discrete chiral symmetry, 

tPix) e^""^^^(x) = (-l)Xx) 

^(x) V^(x)e'""^^ = (-l)>(x), (3) 

where n is an integer. When we introduce a collective field 9{x) via a chiral rotation 

^(x) = e'^(")^^x(^), (4) 
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this degree of freedom is therefore only defined globally modulo nir, i.e, the gauge fixing S- 
function must be periodic. If, in the notation of we take the case A = 0, this gauge-fixing 
function is = 9{x)/'k, and hence 

5{^{x)) = 6{^{x)-n). (5) 

We can give a convenient representation of such a globally periodic functional 5-function by 
a Fourier transform: 



5($(a;)) = J^b]kexp J d^xb{x)^{ 



X] 



(6) 



A:=— oo ' 

where the functional integral is performed over all 6's satisfying the constraint 

1 



TT 



d xb{x) = k, (7) 



where k is an integer. 

The interpolating gauge-fixing function reads in detail, for general A |]l|: 

^x) = aT dCxiOiui.xiO + ^—^oix), (8) 

J-oo TT 

where, for convenience, we have shifted the arbitrary lower limit of the line integral all the 
way to — C)0.[] The periodicity of the 5-function constraint is not immediately obvious if we 
choose A 7^ 0, since under the discrete chiral transformation (j^), the contribution from the 
^-term above now becomes —(1 — A)n. However, the line integral also transforms under the 
discrete chiral rotation, 



A 



r dCxiuX ^ A r dCxiuX - riA, (9) 

J — oo J — oo 



as can be seen by choosing an extremely slowly varying function, a{x), tending to the limit of 
UTT everywhere (except at — oo), and then performing the corresponding regularized discrete 
chiral rotation 

^^^e-W75^ ^ ^ e^<^(^)^^X- (10) 
Under this transformation, the line integral transforms on account of the axial anomaly: 

A r dCx^uX ^ A r dCxiuX--a{x) ^ A f dCxiuX-nA. (11) 

J —oo J —oo TT J — oo 

The two pieces therefore add up, leaving us again with the result (|^). 

Using the Fourier representation (^), we can rewrite the exponent as an interaction term 
coupling the fermions and the boson to a vector potential B^{x) defined by 

B,{x) = e,, [ d'ybiy) f dCS^'\^-x). (12) 

J J —oo 



^In a manner somewhat reminiscent of the prescription used by Mandelstam in relating boson and 
fermion operators. Here, the point — oo should simply be viewed as any point "outside" the space-time region 
we are considering. A more rigorous treatment could either keep xq, the lower limit of the line integral finite, 
and then modify the analysis by the inclusion of zero modes, - or choose appropriate boundary conditions 
in a finite volume. 
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The Pontryagin index density of this vector potential is 



-e'^'d^B^ix) = --bix); (13) 

TT 71 



by the requirement (0), this imphes 



- / (fxt^"'d^B^{x) = k (14) 



7T 



where the integer k is the instanton number of the gauge potential 5^. In general, this gauge 
potential can be decomposed as (see, e.g., ref. |TU|) 



B^{x) = kC^{x) + e^,d''b{x) + d^^{x), (15) 

but we can always remove the third term by a local phase rotation of the fermion fields. 
This simply corresponds to choosing a gauge for B^{x). The field C^{x) is a background 
which is only constrained to have the volume integral of 7T~^e^vd'^C^{x) fixed, i.e. it carries 
topological number 1. We can choose such that d^Cy — dyC^ is constant [|1T]. The relation 
of the measure defined by (§) to those in terms of /c, 6, respectively -B^ is given by 



-j-OO -{-"OO -|-oo 

E jmu= E jv'[h]\det\d')\= E jv'[h]V'[^]5[^)\det\d')l (16) 

K=—oo k=—OD k=—oo 

where the prime denotes that the zero-mode sector (with respect to d"^) is excluded. The 
last expression is nothing but the measure for the vector field 5^ in the decomposition (^5|) 
for the gauge ip = 0. 

Let us now add a term 

^A^r = -^,B,^{x)B^^{x) (17) 

to the action. The field strength is B^i,{x) = d^B,^ — d^B^. Surprisingly, the dimensionful 
coupling constant g will turn out to play the role of an ultraviolet cut-off. The effect of 
adding such a term is to smear the ^-function of the gauge-fixing into a Gaussian. The 
advantage is that, before integrating over 9, our action now contains one part which is the 
Schwinger model coupled to an external source. We can then directly make use of results 



which have been established for that model [T^. Later, we of course send the cut-off g to 



infinity, thereby reducing the gauge-fixing function to the original ^-function constraint. 

The case A = 1: Bosonization 

Putting everything together, we can now represent our gauge-fixed Lagrangian as 



+M(e-''' - lY^^ - — Miie^'' - 1) + M' (e 
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The additional terms depending on M± result from the introduction of the field 6 via a 
chiral transformation. The corresponding Jacobian of such a transformation, calculated 
with a Pauli-Villars regularization, contains also mass-dependent terms (see |]T^, |I[|). With 
the constants Cj, ki obeying the relations 

Clh + 02^2 = , Ci + C2 = 1, (19) 

the cut-off-dependent constant ki(A) is defined as 

K,{A)=AJ2c^h\ogkl (20) 



we use a 



For convenience, and in order to compare our results with the work of Dorn |[13[ 
scheme where Ci = C2 = ^, ki = — ^2 = 1. With this choice ki(A) = 0. 

It should be noted that the terms M^(e^^*^ — 1) in eq. ([T8|) arise as a consequence of this 
particular regularization scheme. As we shall see later, they correspond to contact terms of 
the ^-function kind in certain Green functions. Such (5-function contributions are not naively 
expected to occur in a regularized theory. However, Pauli-Villars regularization is defined 
here as a subtraction scheme of connected n-point functions (see [0), not as a regulariza- 
tion of the free propagator itself. This leads to (5-function contributions when calculating 
quantities involving more than one loop, like the two-loop contribution M^{e^^^^ — 1) to the 
Jacobian mentioned above. It will turn out that these contact terms proportional to 
appear in the bosonized version of the generating functional as well. One could try to avoid 
the occurrence of these terms by using a different regularization method; however, some very 
pleasant features would then have to be given up. In particular, this Pauli-Villars regulariza- 
tion guarantees that independent (commuting) functional derivatives of Z[Vf^, A^j^, M±] yield 
the same result, irrespective of the order in which they are taken, as was shown in ref. |]12| 



The Lagrangian (|T8|) looks almost like a Schwinger model coupled to sources M±{x)e'^'^^^^^\ 
except for the derivative coupling of 6 to the fermions^. This coupling can be removed by 
the following set of transformations: 

9 {l + a)-^9. (21) 

These transformations have to be carried out in the order shown. Here, a is a non-local 
operator: 

a = -d^ (^d' + . (22) 

As a consequence the kinetic term for 6 acquires an additional term, which leads to a Pauli- 
Villars regularized propagator with regulator mass The Lagrangian reads as follows: 

^ = ^Sch + ^int + 



^And of course the additional part of the Lagrangian depending only on 0. Since is a dynamical 
field, the full theory is not completely decoupled into separate sectors. But as we shall henceforth perform 
a perturbative expansion in the terms xM±(a;) exp[±2i0(x)], this has no consequences for the following 
analysis. 
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Ce = +j^d,el^d'' + ^y^e-^i^Ml{e^^'-l) + M^_{e-'^'-l)j. (23) 

Now we are in a position to derive an effective bosonic theory in terms of the field 6. Let 
us therefore treat Cint as a perturbation of the Schwinger model part, i.e. 

(exp^^i J (fxCintix)^) =exp(^i J £x{Cint{x)) + ^— J (fx J £y{Cint{x)Cint{y))c + ■ (24) 

The expectation values have to be taken with respect to x, X aiid B^. The expression (. . .)c 
means the connected part of the correlation function. The resulting effective bosonic action 
reads 

Cm^ = M^{x)e'''^^\xix)P+x{x)) + {+ ^ -) 

+ ^|rf2yM+(x)e^^^(^)(x(x)P+x(a:)x(l/)P+x(l/))cM+(y)e2^^(^) 

+ U d'yM4x)e'^'^^\x{x)P+x{xmy)P-X{y))cM.{y)e-'^'^y^ 



2 

+ (+--) 

+ . . . (25) 

All constants independent of 6 or M± have been dropped. They are unimportant for our 
discussion because they only contribute to the normalization of the functional integral. 

We can now use known results for the Schwinger model, in particular the cluster decom- 
position property [1^. First of all, for the chiral condensates we obtain the well-known 
expression 

{XP±X) = (26) 
y 7r47r 

where 7 is Euler's constant. For the 2-point functions of fermion bilinears we get (in 
Euclidean space, and including the contributions from the Pauli-Villars regulators for the 
fermions [|l^, 0): 

{xix)PM^)x{y)P±xiy))c = j^5ix -y)+ (^-g-)' (^e-'^°(*l^-^l) - l), (27) 

and 

{xix)PM^my)P^xiy))c = (-^)'(e'^°^*l^-^l^e(e^A|x - - 1) - l), (28) 

valid up to terms of 0{A~^). Here Kq is a modified Bessel function of the second kind. Of 
course we take the limit where A is sent to infinity, and as a consequence the step-function 
0(e'^A|a; — y\ — 1) above simply equals unity for physical distance^. The scheme-dependent 
6{x — y) contribution from {x{.x)P±x{.x)x{y)P±x{y))c to the effective action of eq. ( p5|) 
just cancels the (l/87r)M-|-(x)^e^^*^*-^-' -terms of Cg. This cancellation is already one highly 
non-trivial step towards identifying the bosonization relations. 

^That is, distances larger than 1/A. 
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The cluster decomposition property (shown in |]10| to hold even in this massless Schwinger 
model in which we are computing averages) states that the connected part of the correlation 
functions vanishes for distances \x — y\ ^ This can be seen quite easily from the 

above expressions, taking into account the properties of the Bessel function Ko{z) at large 
argument. If we take the limit (7 — > cxo in order to recover the old 5-function gauge, all 
these connected correlation functions vanish for all non-vanishing distances. It may appear 
surprising that the gauge parameter g plays such a crucial role in simplifying the analysis; 
after all, any physical answer is independent of this parameter. However, although we may 
recover the same final results also for finite g, the corresponding equivalent bosonic theory 
will in general be highly non-local, and not very illuminating. As we shall see below, it 
is only in the limit g ^ 00 that we recover the standard local bosonized action. So the 
most convenient bosonization gauge remains the (5-function choice of which corresponds 
precisely to (7 — 00. 

A note should also be made here concerning the distance scales involved at this point. 
We started out with an ultraviolet cut-off A from the Pauli-Villars regularization of the 
original fermion theory. Through the modified gauge-fixing function [the addition of a kinetic 
energy term for the gauge potential (P!7|)], we introduced what turned out to be a Pauli- 
Villars regulator for the boson field 6{x). Finally, we also need an infrared cut-off fi for the 
computation of certain Green functions (see below). Although both ultraviolet cut-offs A and 
g are eventually taken to infinity, we are always performing the analysis in a certain distance 
regime set by these two scales. First of all, from the beginning we clearly restrict the whole 
analysis to distance scales |a; — y| ^ 1/A. The cluster decomposition sets in for distances 
\x ~ y\ > 1/g. Taking the limit g ^ 00 guarantees that all higher connected correlation 
functions vanish at all physical distances. However, if we wish to retain a finite cut-off g, 
then all terms in the expansion contribute to the effective bosonic theory, suppressed only 
by powers of g~^. 

Finally, adding all the above ingredients and rescaling 9 by a/tt, we arrive at a bosonized 
theory described by a Lagrangian 

C,os = ^,d,e (d' + d^e + M+^e^'V^^ + M_ ^e"2.0F^ _ ^(m^ + m!), (29) 

where the limit (7 ^ 00 is to be taken.Q This coincides with the result of Dorn [|1^ if we 
identify g/\/TT with his Pauli-Villars regulator mass Ms- Indeed, the kinetic term for 9 
is precisely the corresponding Pauli-Villars regularized kinetic energy. Taking the sources 
M± = m to be constant, this is the Sine-Gordon action with a Pauli-Villars mass g/ \/t^ and 
one particular normal-ordering prescription in the operator formalism . We hope that the 
present derivation of the bosonization relations (1) for fermionic mass terms has underlined 
some of the subtleties behind the statement that the massive Thirring model is to be viewed 
as equivalent to the 2D Sine-Gordon theory. But the path-integral manipulations of the 
present paper can also be viewed as a much simplified derivation of the classic result of 
Coleman ISll. 



■^It should again be emphasized that there are corrections to this simple effective action, down by powers 
of ■ The shown terms should really be viewed as the result of the limit g — > 00. Everything is here 
expressed in terms of bare quantities. 
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Conclusion 



We have shown how to regain the usual mass-term bosonization identities by means of the 
gauge $ with A = 1. The case A = of course trivially yields the purely fermionic for- 
mulation. As in previous work , the equivalent bosonic Lagrangian has been derived only 
through its perturbative expansion. However, in contrast with ^ and other similar ap- 
proaches using path integral methods 0, we do not explicitly need an order- by-order com- 
parison in the perturbative expansion. Using the established cluster decomposition property 
of the Schwinger model [which a priori would seem to be unrelated to the present treatment 
of the ungauged action (H)], we need only the first two orders of the expansion (|2^) . All 
higher orders vanish. Thus, although we have not yet demonstrated a truly non-perturbative 
boson-fermion equivalence for these mass terms, we feel nevertheless that the present deriva- 
tion represents a substantial improvement. 

As for the restriction at intermediate steps to very definite distance scales, we see here an 
effect which is bound to occur if the same collective field technique is used to derive effective 
Lagrangians (partly bosonized, or not) in higher dimensions. The analysis of was in that 
sense rather particular in that for fermionic couplings to just vector and axial vector currents 
in two dimensions, no true field theoretic regularization was required beyond the calculation 
of functional Jacobians. Here we have seen the more general machinery at work. If we insist 
on a finite cut-off A, the fermionic theory (0) can only be partly be bosonized in a local 
fashion, with rather unpleasant non-local corrections of 0{A~^). In addition, only in the 
limit (7 ^ oo do we immediately achieve a local action in the bosonic representation.^] 

It may be worthwhile to add a few comments on the computation of chiral condensates, 
and in general fermionic correlation functions, of the ong'ma/ fermionic theory (^). Even for 
finite g the physical quantities should not depend on this parameter. Thus, if one calculates 
expectation values (or Green functions) in the effective ^-theory, occurrences of g should 
cancel with contributions from the ^-integration. Indeed, this is what happens. Take as an 
example the chiral condensate {ipP±ip) in a free theory, given by a functional derivative with 
respect to M±, which is subsequently set equal to zero. The integration over 9 has to be 
done with respect to the kinetic term defined in Cq, even for finite g. This yields 

iijP^^P) = -^{e'^') = hm 4£_e-^o(.|.|)+Xo(;^N) ^ (33) 

where /i is an infrared cut-off mass. This means that in going from the formulation in terms 
of rotated fermions x? X to the physical fermion fields and gjs/n is replaced by the 
infrared cut-off /i, which is finally taken to zero. One then recovers the correct result for the 
chiral condensate of a free massless theory. 
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